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Using the numerical solution of the nonlinear Schrodinger equation and a variational method 
it is shown that (3+l)-dimensional spatiotemporal optical solitons, known as light bullets, can be 
stabilized in a layered Kerr medium with sign-changing nonlinearity along the propagation direction. 

PACS numbers: 42.65.Jx, 42.65.Tg 
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I. INTRODUCTION 

After the prediction of self-trapping [i| of an optical 
beam in a nonlinear medium resulting in an optical soli- 
ton 0, y] I there have been many theoretical and exper- 
imental studies to stabilize such a soliton under differ- 
ent conditions of nonlinearity. A bright soliton in (1+1) 
dimension (D) in Kerr medium is unconditionally sta- 
ble for positive or self-focusing (SF) nonlinearity in the 
nonlinear Schrodinger equation (NLS) ^J. However, in 
(2+l)D in homogeneous bulk Kerr medium one cannot 
have a stable soliton-like axisymmetric cylindrical beam 
0,13; H- Also, in (3-|-l)D in such a medium one cannot 
have a stable optical wave packet that remain confined in 
all directions. Such a confined wave packet in (3-|-l)D is 
often called a light bullet and represents the extension of 
a self-trapped optical beam into the temporal domain Q ■ 
If the nonlinearity is negative or self-defocusing (SDF), 
any initially created soliton spreads out in both (2-|-l)D 
and (3-|-l)D ^]. If the nonlinearity is positive or of SF 
type, any initially created soliton is unstable and even- 
tually collapses p\. 

Recently, through a numerical simulation as well as 
a variational calculation based on the NLS it has been 
shown that the axisymmetric cylindrical beam in (2-|-l)D 
can be stabilized in a layered medium if a variable non- 
linearity coefficient is used in different layers [3, ISl- ^ 
weak modulation of the nonlinearity coefficient along the 
propagation direction leads to a reasonable stabilization 
in (2-|-l)D [g. A much better stabilization results if the 
Kerr coefficient is a layered medium is allowed to vary 
between successive SDF and SF type nonlinearities, i.e., 
between positive and negative values j2j. However, it has 
been shown that such a modulation of the nonlinearity 
coefficient in a Kerr medium should fail to achieve stabi- 
lization of a light bullet 9] or a general three-dimensional 
soliton [lOj. 

As the stabilization of a light bullet is of utmost inter- 
est, we revisit this problem and find, to great surprise, 
that a spatiotemporal light bullet can be stabilized in a 
layered Kerr medium with sign-alternating nonlinearity 
along the propagation direction. 

Although, the present work is of interest from a the- 
oretical point of view, it also has phenomenological or 
experimental consequences. Recently, it has been em- 
phasized ill]] that large negative values of the Kerr coef- 
ficient can be created by using the cascading mechanism 
with a large phase-mismatch parameter. It has also been 



suggested T^l that a layered medium with alternating 
sign of nonlinearity can be created with the technique 
of mesoscopic self-organization. Hence a stabilized light 
bullet can be experimentally realized in the future. 

To stabilize a soliton in a SF homogeneous bulk Kerr 
medium, the repulsive kinetic pressure due to the Lapla- 
cian operator in space and time in the NLS has to balance 
the attraction due to nonlinearity. For a light bullet of 
size L, kinetic pressure is proportional to L~^ whereas 
attraction is proportional to L~^ in (D -|- 1) D. The ef- 
fective potential, which is a sum of these two terms, has 
a confining minimum only for _D = 1 leading to a stable 
soliton ;^13]. Using a variational method we find that a 
layered Kerr medium with sign-changing nonlinearity in 
(3+l)D can lead to an effective potential with a mini- 
mum which can stabilize the solitons. 

In Sec. II we present a variational study of the problem 
and in Sec. Ill we present a complete numerical study. 
Finally, in Sec. IV we give the concluding remarks. 
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II. VARIATIONAL CALCULATION 



For anomalous dispersion, the NLS can be written as 



»^ + ^V2+7(z)|u(r,z)|2 



u(r,z) = 0, (1) 



where in (3+l)D the three dimensional vector r has space 
components x and y and time component t, and z is 
the direction of propagation. The Laplacian operator V^ 
acts on the variables x, y, and t. In (2+l)D in Eq. (^ 
the vector r could have components x and y ox x and t, 
while z continues as the direction of propagation. The 
nonlinearity coefficient 7(2) in a layered Kerr medium 
is piecewise continuous and can have successive positive 
(SF) and negative (SDF) values 7+ and 7_ in layers of 
width d. The normalization condition is /(ir|u(r, z)p = 
P, where P is the power of the optical beam ^^,{J\■ 

For a spherically symmetric soliton in (3-|-l)D, 
w(r, z) = C/(r, z). Then the radial part of the NLS (dJ 
becomes [3 
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U{r,z)^0. (2) 



In the following we consider variational and numerical 
solutions of Eq. ((SJ. 



First we consider the variational approach with the 
foUowingtrial Gaussian wave function for the sohition of 
Eq. © |l[l3 



U{r,z) = A^(z)exp 



2R^{z] 



+ —b{z)r + ia{z) 



(3) 



where N{z), R{z), b{z), and a{z) are the soHton's am- 
pUtude, width, chirp, and phase, respectively. In Eq. Q 
in (3+l)D N{z) = P^/'^/[tt^/*R^/'^{z)]. The trial func- 
tion lO satisfies (a) the normalization condition 0, 13 ^^ 
well as the boundary conditions (b) U{r, z) -^ constant 
as r ^ and (c) \U{r, z)\ decays exponentially as r — > oo 



The Lagrangian density for generating Eq. |(21) is [3 



C{U) 



dU 
dz 



U* 



dz 



U 



dU 
dr 



l\U\' 



(4) 



The trial function ||2J) is substituted in the Lagrangian 
density and the effective Lagrangian is calculated by in- 
tegrating the Lagrangian density: L„g — J £{U)dr. The 
Euler-Lagrange equations for R{z), b{z), and a{z) are 
then obtained from the effective Lagrangian in standard 
fashion [3, M, ^M Ehminating a{z), the equations for 
b{z) and R{z) in (3-f 1)D can be written as 



dR{z) 

dz 

db{z) 

dz 



R{z)b{z) 
1 



R\z) 



b\z) 



^{z)P 1 



(5) 
(6) 



From Eqs. ((SJ and © we get the following second-order 
differential equation for the evolution of the width 



d^Rjz) 
dz^ 



l{z)P 



R^{z) 



(7) 



Here we re-visit the stability condition of light bullets 
of Eq. (O for 7(z) = 70 + 7i(z), where 70 = 50 is a pos- 
itive constant of SF type and 71(2:) is a rapidly varying 
part with zero mean value. We take 71(2) = gi sin(a;z), 
as this is a form that we can integrate easily. We break 
R{z) into a slowly varying part A{z) and a rapidly vary- 
ing part B{z) by R{z) = A{z) + B{z). Substituting this 
into Eq. ^ and retaining terms of the order of lo^"^ 
in B{z) we obtain the following equations of motion for 
B{z) and A{z): 



d^Bjz) 
dz^ 



giPsin(wz) 

2V2^ A\z)' 



(8) 



d^Ajz) 
dzJ 



9oP 



A^iz) 2V2^A'^{z) 



2giP{B{z)sm{ujz)) 
V2^A^z) 



where ( ) denotes time average over rapid oscillation. 
Using the solution B{z) = giPsm{u;z)/[2V27r^uj'^A^{z)], 



the equation of motion for A{z) becomes 

d^Ajz) ^ 1 goP , gfP^ 

dz^ 




The quantity in the square bracket in Eq. Hll|) is the 
effective potential U{A) of the equation of motion 

^(^) = 2l^ - ^7=r^ + 32^^218- (12) 

Stabilization is possible when there is a minimum in 
this effective potential [13 . Unfortunately, this condition 
does not lead to a simple analytical solution. However, 
straightforward numerical study reveals that this effec- 
tive potential has a minimum for a positive goP corre- 
sponding to attraction (SF nonlinearity) with ggP above 
a critical value. For a numerical calculation the quantity 
gP is taken to be of the form gP = g^P + giP sin(a;2:) = 
goP[l + 4sin(f07rz)], so that gi = Ago and ui = IOtt. The 
numerical values for gi and uj are taken as examples, oth- 
erwise they do not have great consequence on the result 
so long as lu is large corresponding to rapid oscillation. 
In Fig. 1 (a) we plot the effective potential U{A) vs. 
A for goP = 30, 100, 200, and 500 for w = IOtt. For 
goP = 30 there is no minimum in U{A), whereas a mini- 
mum has appeared for ggP = 100 which becomes deeper 
for goP = 200 and 500. A careful examination reveals 
that the threshold for the minimum in the present case 
is given by goP « 40. Hence in the present case stabi- 
lization is not possible for goP = 30, and it is possible 
for goP > 40. There is no upper limit for 170^ and sta- 
bilization seems possible for an arbitrarily large goP. As 
goP increases the depth of the effective potential in Fig. 
1 (a) increases and consequently, it is easier to stabilize 
a soliton. As the first and the third terms on the right 
hand side (rhs) of Eq. H10|) are positive, no stabilization 
is possible for a negative goP corresponding to repulsion 
(SDF). 

In order to see the effect of the frequency cu on stabi- 
lization, we plot in Fig. 1 (b) the same effective potentials 
of Fig. 1 (a) for uj = 30tt. With the increase of uj the 
effective potentials have become deeper and hence the 
stabilization easier. For goP = 30 we have a minimum 
in Fig. 1 (b) whereas there is no minimum in Fig. 1 
(a). With the increase of u) the threshold value of goP 
for obtaining a minimum has reduced. 

The sinusoidal variations of Kerr nonlinearity as con- 
sidered above for variational study only simplifies the 
algebra and is by no means necessary for stabilization 
of solitons. In the following numerical study we estab- 
lish that a rapid oscillation of the nonlinearity coefficient 
between positive and negative values also stabilizes the 
soliton in (3+l)D. 



III. NUMERICAL CALCULATION 

We solve Eq. |(2Jl numerically using the split-step time- 
iteration method employing the Crank-Nicholson dis- 
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FIG. 1: The effective potential U(A) of Eq. (12) vs. A in 
arbitrary units for goP = 30, 100, 200, 500 (from shallower to 
deeper minimum), for (a) lu = IOtt and (b) lu — SOtt. 



cretization scheme [T^. The time iteration is started 
with the known solution of some auxiliary equation with 
zero nonlinearity. The auxiliary equations with known 
Gaussian solution are obtained by adding a harmonic os- 
cillator potential r^ to Eq. (0). Then in the course of 
time iteration the power P and a positive constant SF 
Kerr nonlinearity 7(2:) = go = 1 is switched on slowly and 
the harmonic trap is also switched off slowly. If the non- 
linearity is increased rapidly the system collapses. The 
tendency to collapse or expand must be avoided to obtain 
a stabilized soliton. 

After switching off the harmonic trap in Eq. ^ and af- 
ter slowly introducing the final power P and the constant 
nonlinearity "f{z) = go = 1, an oscillating Kerr nonlin- 
earity 71 (z) oscillating between zLgi {\gi\ > 1) like a step 
function as z changes by d is introduced on top of the 
constant nonlinearity. The overall Kerr nonlinearity now 
has successive positive and negative values 7+ = 1 + gi 
and 7_ = 1 — gi of equidistant layers of width d in z 
direction. A stabilization of the final solution could be 
obtained for a suitably chosen 7± and a small d. If the SF 
power after switching off the harmonic trap is large com- 




FIG. 2: Wave function \U{r)\ of Eq. (2) of stabilized light 
bullet of power P — 682.6 and (7(2:)) = go = 1 and 71(0) 
oscillating between 7+ = 5 and 7_ = —3 in successive layers 
of width d = 0.1 in the z direction. The wave functions are 
shown at positions z = 0, 50, 100, 150, 200, 250, 300, 350 and 
400. 



pared to the spatiotemporal size of the beam the system 
becomes highly attractive in the final stage and it even- 
tually collapses. If the final power after switching off the 
harmonic trap is small for its size the system becomes 
weakly attractive in the final stage and it expands. The 
final nonlinearity has to have an appropriate intermedi- 
ate value, decided by trial, for final stabilization. The 
stabilization could be obtained for a large range of val- 
ues of 7± and d. After some experimentation with Eq. 
(O we opted for the choice 7+ = 5, 7_ = —3, and d — 0.1 
in all our calculations. 

In the present scheme of stabilization the nonlinear- 
ity rapidly fiuctuates between appropriate positive (SF) 
and negative (SDF) values in the propagation direction. 
For positive nonlinearity, the system tends to collapse, 
whereas in the SDF regime it tends to expand to infinity. 
If the nonlinerities are appropriate, the collapse in the 
first interval is exactly compensated for by the expan- 
sion in the next interval and a stabilization of the system 
is obtained. Obviously, the system will be more stable 
when the intervals are small so that the fiuctuation of the 
system around a stable mean position is small. Conse- 
quently, the system remains virtually static and the very 
small oscillations arising from collapse and expansion re- 
main unperceptable. 

Although, for the sake of convenience we applied a har- 
monic trap in the beginning of our simulation, which is 
removed later with the increase of nonlinearity, this re- 
striction is byno means necessary for stabilizing a soliton. 
Saito et al. 13] used a qualitatively similar, but quantita- 
tively different, procedure for stabilization in the context 
of Bose- Einstein condensation in two dimensions. The 
procedure of Saito et al. could also be applied sucess- 
fuUy in the present context. 

Now we turn to a numerical investigation of Eq. 1^)1. 
The results are shown in Fig. 2 for the (3-t-l) dimensional 



soliton, where we plot the radial part of the wave function 
for different z for a power P = 682.6. A finetuning of the 
power was needed for the stabilization reported in Fig. 
2. 

In Fig. 2 the narrow spread of the wave function over 
the large interval of z shows the quality of stabilization. 
The results at intermediate z lie in the region covered 
by the plots. The plot of the full wave function at dif- 
ferent z on the same graph clearly shows the degree of 
stabilization achieved. The stabilization seems to be per- 
fect and can easily be continued for longer intervals of 
z by increasing the power. In Ref. layers of width 
d — 0.001 was employed for stabilization in (2+l)D. The 
present stabilization is obtained with a much larger width 
d = 0.1, which make the present proposal more attractive 
from a phenomenological point of view. The stabilization 
can only be obtained for beams with power larger than a 
critical value. Numerically, we found it was easier to ob- 
tain stabilization of beams with power much larger than 
the critical value. In (3-)-l)D good stabilization could be 
obtained for much larger power: the power employed in 
Fig. 2 was 682.6, whereas the critical variational power 
for stabilization obtained in Fig. 1 is about 40. 

Using a variational procedure alone, not quite identi- 
cal with the present approach, in the context of Bosc- 
Einstein condensation AbduUaev et al. [3 also had 
found that a stabilization of a soliton could be possi- 
ble in (3-|-l)D via a temporal modulation of the non- 
linearity. However, they confirmed after further analyt- 
ical and numerical study that such a stabilization does 
not take place in (3-|-l)D. Saito et al. 0| and Towers 
et al. 13, on the other hand, are silent about the pos- 
sibility of the stabilization of a soliton in (3-|-l)D. We 
point out one possible reason for the negative result ob- 
tained by AbduUaev et al. 9] in (3-f 1)D. The nonlin- 
earity parameter XqN = ^/2'k^^'^A with A = 1 used in 
Ref. PJ for stabilizing a soliton in (3+l)D is much too 
small (smaller than the threshold discussed in Sec. II). 
Comparing Eq. (2) of Q with our Eq. ^ we find that 
the above value of nonlinearity corresponds in our no- 
tation to goP = XoN/2 = 7r3/2/V2 « 3.9, whereas the 



present variational threshold for obtaining a stabilized 
soliton isQoP ~ 40. The very small nonlinearity used 
in Ref. U is most possibly the reason for the negative 
result obtained there. More recently, Montesinos et al. 
[iQ] have also confirmed the conclusion of Ref. 9] that 
no stable three-dimensional soliton could be obtained by 
a variation of Kerr nonlinearity. However, they did not 
give details of their study, which led to this conclusion, 
for a comparison. 



IV. CONCLUSION 

In conclusion, after a variational and numerical study 
of the NLS we find that it is possible to stabilize a spa- 
tiotemporal light bullet in (3+l)D by employing a layered 
Kerr medium with a sign-changing nonlinearity along the 
propagation direction. From a variational calculation we 
show that a oscillating Kerr nonlinearity produces a min- 
imum in the effective potential, thus producing a poten- 
tial well in which the soliton can be trapped. The present 
stabilized soliton is a slowly collapsing Townes soliton 
[1] with large power. The oscillating Kerr nonlinearity 
stops the collapse and enhances the lifetime of the soli- 
ton greatly. This is of interest to investigate if such light 
bullets could be created experimentally. 

Apart from optics, such stable three-dimensional soli- 
tons can be realized experimentally in Bose-Einstein con- 
densates (BEC), where a Feshbach resonance could be 
used to generate an oscillating nonlinearity or an oscillat- 
ing effective interatomic interaction via the modulation 
of an external background magnetic field M, Il3l ■ The 
stabilization of such BEC solitons in two |9lll0l|13l| and 
three [I^ dimensions is already under investigation. 
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